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ABSTRACT 

Multi-band Hubble Space Telescope photometry reveals that the main se- 
quence, sub-giant, and the red giant branch of the globular cluster NGC 6752 
splits into three main components in close analogy with the three distinct seg- 
ments along its horizontal branch stars. These triple sequences are consistent 
with three stellar groups: a stellar population with a chemical composition sim- 
ilar to field halo stars (population a), a population (c) with enhanced sodium 
and nitrogen, depleted carbon and oxygen and enhanced helium abundance 
(AY ~0.03), and a population (b) with an intermediate (between population 
a and c) chemical composition and slightly helium enhanced (AY ~0.01). These 
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components contain ~25% (population a), ~45% (population b), and ~30% 
(population c) of the stars. No radial gradient for the relative numbers of the 
three populations has been identified out to about 2.5 half mass radii. 

Subject headings: globular clusters: individual (NGC6752) - - stars: Popula- 
tion II 



1. Introduction 

It is now widely accepted that many (if not all) globular clusters (GCs) are made up 
of multiple populations of stars. Since the 70's, it was well known that stars within the 
same cluster have light-element abundance variations (e.g. Kraft 1979, Norris & Freeman 
1979), but, at that time it was not clear whether these 'abundance anomalies' were due 
to internal mixing or differences in the primordial composition or a combination of these 
effects. More recently, high-resolution spectroscopy has revealed the presence of well-defined 
correlations among some light-elements abundances (e.g. Kraft et al. 1992, Sneden et al. 
1994, Ramirez & Cohen 2002), including the anticorrelations between Na and O, and Mg 
and Al, which indicate that material has been processed via high-temperature proton capture 
nucleosynthesis (Denisenkov & Denisenkova 1989). 

The fact that the same light-elements variations have also been observed in unevolved 
cluster stars (e.g. Briley et al. 1994, 1997, Cannon et al. 1998, Gratton et al. 2001), whose 
internal temperatures do not allow high-T proton captures, and in fully-convective low- 
mass M-dwarfs (Milone et al. 2012) suggests that these stars were born with these chemical 
peculiarities imprinted in the matter from which they formed (Cottrell & Da Costa 1981, 
see Gratton et al. 2004 for a review). 

High-precision Hubble Space Telescope (HST) and ground-based photometry has shown 
that several GCs host multiple main sequences (MSs) including ui Centauri, NGC2808, 47 
Tuc, NGC 6752, and NGC 6397 (Anderson 1997, Bedin et al. 2004, Piotto et al. 2007, Milone 
et al. 2010, 2012a,b), which have been associated with stellar populations with different 
helium abundances (D'Antona & Caloi 2004, Norris 2004, Bedin et al. 2004, Piotto et al. 
2005, D'Antona et al. 2005). Multiple stellar populations have also been detected in the 
color-magnitude diagram (CMD) from the presence of multiple sub-giant branches (SGBs, 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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Milone et al. 2008, Anderson et al. 2009, Marino et al. 2012, Piotto et al. 2012) or multiple or 
spread red-giant branches (RGBs, e.g. Grundahl et al. 1998, 2000, Yong et al. 2008, Marino 
et al. 2008, Lee et al. 2009). 

Stellar evolution models predict that high-temperature H-burning through the CNO 
cycle should result in an increase of the N abundance, at the expenses of C and O, and of an 
increase in the helium fraction. Multiple stellar populations with different helium content 
could also account for the HB morphology of some GCs in which the bluer HB sequences can 
be associated with the presence of He-rich stars (e.g. D'Antona et al. 2002, D'Antona & Caloi 
2008, Busso et al. 2007, Cassisi et al. 2009, Catelan, Valcarce & Sweigart 2010, D'Alessandro 
et al. 2011). Indeed, clear evidence of the connection between the HB morphology with the 
multiple populations comes from Marino et al. (2011) who have found that blue-HB stars of 
the GC M4 are all Na-rich and O-poor (hence He-rich), whereas red-HB stars are primarily 
Na-poor and O-rich (He-poor) (see also Norris 1981, Smith & Norris 1993). Similar results 
have been found in NGC2808 (Gratton et al. 2011). 

This work adds yet another cluster (NGC6752) to the growing list of clusters with 
photometric and spectroscopic evidence of multiple sequences along the RGB, MS, SGB, and 
HB. We use HST filters covering a wide range of wavelengths to study the multiple stellar 
populations in the GC NGC6752. The presence of star-to-star light-element variations in 
the cluster has been widely reported in the literature (Norris et al. 1981, Grundahl et al. 
2002, Yong et al. 2003, 2005, 2008, Carretta et al. 2007). 

Photometric evidence for three populations of stars along the RGB of NGC 6752 with 
different Mg, Al, Mg, Si, Na, and O content was early identified by Grundahl et al. (2002) 
and Yong et al. 2008 (see also Carretta et al. 2012). These authors found that the Stromgren 
photometric index c\ correlates with nitrogen abundance in stars both brighter and fainter 
than the RGB bump, and suggest that the observed photometric scatter is due to stellar 
populations with different N abundance (see also Milone et al. 2010, Kravtsov et al. 2011, 
Sbordone et al. 2011). 

The HB of NGC 6752 revealed also a complex structure, with two discontinuities that 
define three HB segments (Momany et al. 2002, 2004). Villanova et al. (2009) analyzed 
spectra of seven HB stars with effective temperature ~8000< T cff <9000 K and found that 
six of them have a chemical composition similar to field-halo stars, including helium. A 
recent photometric analysis of data collected with the Advanced Camera for Survey (ACS) 
on board HST also showed that NGC 6752 has a broadened MS with some indication of a 
MS split (Milone et al. 2010) thus suggesting that its stellar populations also have different 
helium content. 
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The paper is organized as follows. In Sect. [2] we present the data and describe the data 
reduction. In Sect. El HI and[5]we study the triple MS, SGB, and RGB, respectively. An effort 
was made to disentangle the multiple-populations in each of these three evolutionary sequence 
separately. In Sect. Owe explore possible theoretical interpretations and estimate the helium 
difference among the three stellar populations. The study of the radial distributions of the 
various stellar populations is also undertaken. Finally, a summary and some discussion follow 
in Sect. 

2. Observations and data reduction 

This work makes use of two data sets. For the central regions of the cluster we used both 
archival and proprietary material collected with the UV- Visual (UVIS) and infra-red (IR) 
channels of the Wide Field Camera 3 (WFC3), and the wide-field channel (WFC) of the Ad- 
vanced Camera for Surveys (ACS) mounted at the HST . Proprietary images (under program 
GO-12311, PI: Piotto) were collected in two 1-orbit visits taken at two different orientations 
the first one on March 23, and the second one on April 3, 2011, and consist of 12x360 s 
images in camera/channel/filter WFC3/UVIS/F275W, and 2x50 s in WFC3/UVIS/F814W. 
We also used the photometric catalogs presented by Anderson et al. (2008) and obtained 
from ACS/WFC images taken under GO-10775 (PI: Sarajedini, see Sarajedini et al. 2007). 
The archive HST material is described in table [2] and consists of images taken through four- 
teen filters spanning a wide spectral range, from the ultraviolet (F225W) to the infrared 
(F160W). Footprints of HST images are shown in Fig. [TJ 

Photometry and relative positions of stars in HST /WFC3 images were extracted with 
the software tools described in Bellini et al. (2011), mostly based on the software described in 
Anderson et al. (2006). The photometry was calibrated onto the Vega-mag system following 
the procedures given in Bedin et al. (2005), and using encircled energy and zero points 
given at STScI's web pages. Star positions were corrected for geometric distortion using the 
solutions given by Bellini & Bedin (2009) and Bellini et al. (2011) for WFC3/UVIS, Anderson 
& King (2006) for ACS/WFC, and Anderson et al. (in preparation) for WFC3/NIR. 

To study the external regions of the cluster, we made use of the ground-based photomet- 
ric catalog published by Grundahl et al. (2002). They have been obtained with the 1.54 m 
Danish telescope at La Silla (Chile) through the Stromgren filters u, v, b, y, and cover a field 
of view of ~6x6 arcmin centered on the cluster. They were reduced following the method 
outlined in some detail by Stetson (2005). 

The stellar catalogs were purged of poorly measured objects using quality indices that 
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our software produces following the procedure that it is described in Milone et al. (2009). 
Finally, the photometry was corrected for zero-point spatial variations following the recipes 
in Milone et al. (2012). 

ARA [arcmin] 

-a -l o l 2 
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O 

-60:00:00 



i RA 

O 

Fig. 1. — Footprint of the HS T fields studied in this paper. The footprint of ACS/WFC, 
WFC3/UVIS, and WFC3/NIR images are colored green, gray, and red, respectively. 



INSTR 


DATE 




NxEXPTIME 


FILTER 


PROGRAM 


PI 


WFC3/UVIS 


Aug 21 2010 




6 x 120s 


F225W 


11904 


Kalirai 


WFC3/UVIS 


Mar 23, Apr 3 2011 




12x360 


F275W 


12311 


Piotto 


WFC3/UVIS 


May 05 2010 




30s + 2 x 500s 


F336W 


11729 


Holtzman 


WFC3/UVIS 


May 05 2010 


2x2s 


+ 2x348s+ + 2 x 880s 


F390W 


11664 


Brown 


WFC3/UVIS 


May 05 2010 




50s + 2 x 700s 


F390M 


11729 


Holtzman 


WFC3/UVIS 


May 05 2010 




90s + 2 x 1050s 


F395N 


11729 


Holtzman 


WFC3/UVIS 


May 05 2010 




40s + 2 x 800s 


F410M 


11729 


Holtzman 


WFC3/UVIS 


May 05 2010 




40s + 2 x 400s 


F467M 


11729 


Holtzman 


WFC3/UVIS 


Aug 7, 21 2010 




12 x 670s 


F502N 


11904 


Kalirai 


WFC3/UVIS 


May 05 2010 




5s + 40s + 400s 


F547M 


11729 


Holtzman 


WFC3/UVIS 


Jul 31, Aug 7, 21 2010 




15 x 550s 


F555W 


11904 


Kalirai 


WFC3/UVIS 


May 01 2010 




30s + 2 x 665s 


F555W 


11664 


Brown 


WFC3/UVIS 


Jul 31, Aug 7, 21 2010 




15x550s 


F814W 


11904 


Kalirai 


WFC3/UVIS 


May 01 2010 




30s + 2 x 495s 


F814W 


11664 


Brown 


WFC3/UVIS 


Mar 23, Apr 3 2011 




2x50 


F814W 


12311 


Piotto 


WFC3/NIR 


May 01 2010 


3x4s + 


3 x49s + 299s+ 2 x 399s 


F110W 


11664 


Brown 


WFC3/NIR 


May 01 2010 


3x4s + 


3 x49s + 299s+ 2 x 399s 


F160W 


11664 


Brown 


ACS/WFC 


May 24 2006 




1x2s + 4x35s 


F606W 


10775 


Sarajcdini 


ACS/WFC 


May 24 2006 




1x2s + 4x40s 


F814W 


10775 


Sarajedini 



NGC 6752 




Table 1: Description of the HST data sets used in this paper. 
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3. A triple main sequence in NGC 6752 




Fig. 2.— m F336W vs. m F336W - m F395N (left panel) and m F275W vs. m F275W - m F336W CMDs 
of NGC 6752. The insets show the Hess diagram of a zoomed section of the MS. 



A visual inspection of the large number of CMDs that we obtain from our data set 
confirms that multiple sequences along the MS and the RGB can be easily identified by 
using different combinations of the F275W, F336W, and F395N filters. The left-hand panel of 
Fig. [2] shows the m F336W versus m F336W — m F3 g 5 N CMD after that the photometric corrections 
and the quality selection described in the previous Section were applied. We note a bimodal 
RGB, and a spread MS. The m F2 7sw versus m F275 w — ^F336W CMD plotted in the right 
panel reveals an even larger number of features, with a possible triple RGB, and a clear split 
MS composed of two distinct components. A narrow red MS, containing about one fourth 
of MS stars, and a more dispersed blue MS. 

In recent papers, we have shown that two-color diagrams obtained by combining a far- 
UV filter (such as F225W or F275W), a near-UV filter (such as F336W) and a visible filter 
(such as F438W) are powerful way for identifying populations of stars with different helium 
and light-element abundances (see Milone et al. 2012a,b for results on 47 Tuc and NGC 6397). 

Motivated by these results, in Fig. [3^ we plot m F275 w — m F336W against m F336 w — ^f4iom 
for MS stars with 19.65 < m F2 7sw < 23.25. Panel b shows the same two-color diagram for 
the RGB stars. 

We also defined the color index c F2 75W,F336W,F4iom=(^F275W - ™F336w) - (™F336W - 
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^F4iom)- Quite interestingly, the m F 275w versus cf275w,F336W,F4iom color-index-Magnitude 
Diagram (CMD) of Fig. [3fc allows us to identify multiple sequences along the entire CMD, 
from the MS to the RGB tip. There is a clear color spread, with the presence of three RGBs 
and two distinct MSs, in close analogy with what was observed in 47Tuc and NGC6397. 
Also the SGB is not consistent with a simple stellar population. In the following we will refer 
to the less-populated MS located on the bottom-left of the two-color diagram, in Fig. 2a, as 
MSa. The MS is analyzed in this Section, while Sects. @] and [5] are dedicated to the SGB 
and the RGB. 




Fig. 3. — mp275w — m F336W versus mF336w — ?^F4iom two color diagram for MS stars with 
19.65 < mp275w < 23.25 (panel (a)), and for RGB stars with mF27sw < 18.0 (panel (b)). 
Panel (c): The m F 27sw versus cf275w,F336W,F4iom diagram for all NGC6752 stars in our 
sample. 

In Milone et al. (2010), we used high-precision photometry of ACS/ HST images to 
search for signs of multiple populations in NGC6752. We found a broadened MS, and 
demonstrated that this broadening is intrinsic. We also noted a possible MS bi-modality, 
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and suggested that the MS split could be due to two stellar populations with almost the same 
age and iron abundance, but different helium content. As demonstrated in the following, 
the data set presented in this paper allows us to identify multiple populations with different 
helium content with an higher accuracy than was possible for Milone et al. (2010), with the 
data available at that time. 

Both theoretical arguments and observations indicate that CMDs with wide color base- 
lines can be very sensitive to helium differences among stars (e.g. D'Antona et al. 2005, 
Piotto et al. 2007). In this context, the m F8 i4w vs - m F275W — m F8i4W CMD shown in Fig. @] 
clearly reveals a bimodal MS. The MS bi- modality is even more evident in the Hess diagram 
plotted in the inset. The two MSs merge close to the turn off, and the MS separation in- 
creases towards fainter magnitudes, from about 0.1 mag at mp8i4w~19.05 up to 0.25 mag 
at m F814W ~20.15. Hereafter, we will refer to the bluest MS of Fig. @]as MSc. We will also 
demonstrate that the MSa and the MSc correspond to different stellar populations. Note 
the different morphology (i.e. different distribution in color of the stars) of the MSs plotted 
in Fig. H] and Fig. [2J right panel. 




F276W JJJ Fai4Vf 

Fig. 4. — m F814W vs. m F2 75w — m F8uw CMD for MS stars of NGC6752. The Hess diagram 
in the inset is a zoom of the region where the MS bimodality is more evident. 

In order to compare the sequences identified in this 'regular m F 8i4W vs. m F27 5w— wifsuw 
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CMD' with those identified in the mF27sw vs. cf275W,F336W,F4iom CMD of Fig. 2c, we color- 
code the MSa stars green, and in blue the MSc stars (see panels (a) and (b) of Fig. [5]). Panel 
(c) of Fig. [5] shows that MSa stars are distinct from MSc, and that MSa+MSc stars are not 
all MS stars in NGC 6752. In the bottom panels of Fig. we identify the MSa and MSc 
stars in both panels, and identify as 'MSb' stars those that are neither MSa or MSc. These 
stars are colored magenta. 

This approach allowed us to demonstrate that the MSa and the MSc correspond to 
two distinct stellar populations, and that NGC 6752 hosts a at least a third MS population, 
labeled as MSb. 




Fig. 5. — Visualization of the procedure to define the sample of MSa, MSb, and MSc stars. 
The red dashed-dotted lines in panels (a) and (b) are used to select MSc and MSa stars that 
we colored blue and green respectively. In the panels (c) and (d) we adopted a magenta 
color to identify stars that belong neither to the MSa nor to the MSc, and hence are part of 
the third MS component (MSb). 
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The method used to estimate the fraction of stars in the MSa is illustrated in Fig. [6j The 
left-hand panel of the figure is a reproduction of the CMD of Fig. [3b with the fiducial line of 
the most populated MS superimposed. The verticalized MS is plotted in the middle panel 
and the right-hand panels show the histograms of the A(cf275w,F336W,F4iom) color distribution 
in five mF275W intervals. We fitted each histogram with the sum of two Gaussians, colored 
green and black. Hereafter, the green color code will be used to highlight MSa stars. From 
the area under the Gaussians we estimate that 25±2% of stars belongs to MSa. The errors 
were computed as the rms of the values obtained for the five bins, and then divided by two 
(i.e. the square-root of the number of bins minus one). 




Fig. 6. — Left panel: Reproduction of the diagram of Fig. [3b, the red line is the fiducial 
of the most populous MS. Middle panel: The same diagram, after subtraction of the color 
of the fiducial line. Right panel: The Acf275w,F336W,F4iom color distribution in five mp275w 
intervals. The gray lines represent the least-square best fits of two Gaussians to the observed 
distribution. 
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The procedure used to obtain the fiducial line and to verticalize the MS has been adopted 
in several previous papers from our group (e.g. Piotto et al. 2007), and will be used several 
times in the remaining part of this paper. Briefly, we adopt the most populous (in this CMD, 
the bluest) MS as reference sequence and draw by hand a first approximation ridge line. We 
also select a color range around this line to include most of the stars on the blue MS. Then 
we divide the reference sequence in 0.15 mag intervals, and calculate the median colors of 
the stars within each interval. These median points are then interpolated with a spline. We 
calculate the spline at the magnitude level of each star and subtracted it from each star's 
color to estimate the MSRL residual for each. We then determined a sigma-clipped mean for 
each magnitude interval, and repeated the procedure several times. The result is a fiducial 
line plotted in Fig. |6j Finally we subtract from the cf275w,F336W,f4iom color of each star the 
corresponding color on the fiducial line at the same mp275w- 

In order to calculate the fraction of MSc stars we followed a recipe similar to the one 
already applied to the MSa, and illustrated in Fig. [71 In the left-hand panel we plotted the 
^F8i4W vs. ^F275W — m F8i4W CMD, while, in the central panel, we show the verticalized CMD, 
obtained with the same procedure as explained above. The color distribution histograms for 
five magnitude intervals are shown in the right-hand panels and are fitted with two Gaussians 
colored blue and black. In the following, the blue color will be used to indicate MSc stars. 
From the area under the best fit Gaussians, we estimate that 31±3% of the total number of 
MS stars belong to MSc. Errors are calculated as described above for the case of the MSa. 
Since we have already estimated that the MSa and the MSc contain 25±2% and 31±3% of 
the total number of MS stars, we can conclude that the MSb is made up of the remaining 
44±4% of MS stars. 

With these identifications and the large number of filters through which we have obser- 
vations, we can analyze the relative location of the three MSs in a large number of CMDs. 
As an example, left-panel of Fig. [H] shows the mp336W versus mp336w — m F39ow CMD for MS 
stars. In the right-panel we plotted the same CMD, but using the same color codes previ- 
ously defined in order to highlight MSa, MSb, and MSc stars. Contrary to what is observed 
in the mp275w — m F8i4W color in Fig. 3 (or in Fig. 4), in this CMD the MSa is bluer than 
the bulk of MS stars, with the MSb and the MSc being almost coincident. 

In order to follow the behaviour of the three MSs in all the possible combinations of 
the photometric bands available in our data set, we followed the approach illustrated in 
Fig. [9j In each panel we show the fiducial lines of the three MSs in the mpsuw versus 
m x - m F814W (or m F814W - m x ) plane, where X=F225W, F275W, F336W, F390M, F390W, 
F395N, F410M, F467M, F502N, F547M, F555W, F606W, F110W, and F160W. For the two 
IR filters, F110W and F160W, we adopted the ttifsmw — color baseline, being in these 
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N 




F275W m F8i4W A(m F275w -m F814W ) A(m F275W -m FB14W ) 



Fig. 7. — m F814W vs. m F27 5w— "^F8i4W CMD and verticalized m F814W vs. A(m F2 75w— m F8i4w) 
diagram (left and middle panels). Right panels show the histogram of the color distribution 
in five m F8 i4w intervals. Gray lines represent the least-square best fits of two Gaussians to 
the observed distribution. 

cases X the redder filter. 

The MSs color properties in the various CMDs can be summarized as follows: 

i) The MSb is typically bluer than the MSa in all CMDs of Fig [9j with the exception 
of the m F814W vs. wi F336W — m F814W and the m F814W vs. m F3 goM — ^F8i4W CMDs, where the 
two sequences invert their relative colors. 



ii) The MSc is bluer than the MSa in all CMDs. 
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Fig. 8. — Left panel: mp336W versus m F 336w — ^F39ow CMD for MS stars. Right panel 
shows an example of the definition of the fiducial lines. The sample of MSa, MSb, and MSc 
stars defined in Fig. [5] are colored green with black shadow, magenta with white shadow, 
and blue with black shadow, respectively. In this case, MSb and MSc fiducials are largely 
overlapping. The same color codes adopted in the previous figures are also used to represent 
the corresponding fiducials. 



iii) The color distance between the MSa and both the MSb and the MSc increases for 
larger color baselines, with the exception of the mp336w — ^F8i4W> the mp390M — m F8i4W, 
and the mp39ow — ^F8i4W colors. 

iv) The color separation between the MSc and the MSb increases with the size of the 
color baseline for all colors studied in this paper. 

In the following, we will use these data to gather information on the chemical composi- 
tion of the three MSs. 



4. Multiple stellar populations along the Sub Giant Branch 



The first photometric evidence of multiple stellar populations along the SGB of NGC 6752 
comes from the recent paper by Kravtsov et al. (2011). Using wide-field ground-based pho- 
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Fig. 9. — MS fiducials in fourteen m F 8i4w versus mx — ttifshw ( or m F8uw ~ m x) CMDs 
(X=F225W, F275W, F336W, F390M, F390W, F395N, F410M, F467M, F502N, F547M, 
F555W, F606W, F110W, and F160W). At the top of each panel we give the color distance 
from the MSa of the other two MSs, measured at 777^^=18.5. The positions of MSa, MSb, 
and MSc at 777^^=18.5 are represented with green, magenta, and blue circles, respectively 
in the inset of each CMD. 

tometry, Kravtsov et al. (2011) identified a spread of ~0.3 magnitudes in U band, with the 
faintest SGB more centrally concentrated than its brighter counterpart. 
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Our multi-color set of CMDs reveals an even more complex picture for the SGB of 
NGC 6752. A visual inspection at the m F336 w versus m F336W — m F814W CMD and the mp275w 
versus cf275w,F336W,F4iom diagram in the upper panel of Fig. [TU1 immediately reveals that, in 
the WFC3/UVIS field of view, there is no evidence for a wide magnitude spread along the 
SGB in F336W (which is the HST analog of the standard U), even though the SGB is clearly 
not consistent with a single stellar population. By analogy with what we did for the MS 
stars, in the lower panels of Fig. [TOj we plot the mp275w — m F336W versus m F336W — m F410 M 
two-color diagram. This diagram shows a multi-modal distribution of SGB stars (see also 
the histrograms of panels f). As shown in the figure, we selected by eye three groups of 
SGB stars that we named SGBa, SGBb, and SGBc and colored green, magenta, and blue, 
respectively. 




Fig. 10. — ^F336W versus m F336W — m F8uw CMD (panel a) and m F 27sw versus 
c F275W,F336W,F4iom diagram (panel b) zoomed around the SGB. SGB stars are highlighted in 
thick black. Panel (c): m F2 7sw — m F336W versus m F 336w — m F4ioM two-color diagrams for the 
stars shown in the upper panels. The red-dashed line is a fiducial line drawn by hand through 
the middle SGB. In the panel (d), three groups of SGBa, SGBb, and SGBc stars are defined 
and color-coded in green, magenta, and blue, respectively. The verticalized m F 275w — ^F336W 
versus A(m F 336w — ^f4iom) is plotted in panel (e), while the histogram of the distribution in 
A 

(^F336W — ^f4iom) is shown in panels (f) for the three quoted m F2 75W ~ m F336W intervals. 
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In order to better understand the properties of these three sequences, Fig. [TT] gives a 4 
x 4 array of CMDs, where stars of the three sequences selected in Fig. [10] are plotted with 
their color code. Figure [TT] shows some significant features of the SGB of NGC6752: 

(i) SGBa stars share some similarities with MSa: they are on average redder than the 
bulk of SGB stars in mp275W-' m F336W, but they become bluer than the remaining SGB stars 
in the other CMDs of the first row of Fig. [TTJ 

(ii) In the F336W band, SGBa stars are typically brighter than the other SGB stars. 
This fact also explains why they appear redder than SGBb and SGBc stars in mp275w — 
^F336W ; and bluer than them in the other CMDs of the third row of Fig. [TTJ 

(iii) In all the CMDs, the SGBb sequence seems to be located between the SGBa and 
the SGBc. 

The analogy in color distribution of the three MSs and SGBs justifies the names that we 
gave to these sequences, which explicitly want to suggest that the SGBa, SGBb, and SGBc 
represent the continuation along the CMD of the MSa, MSb, and MSc, respectively. 

5. Multiple stellar populations along the Red Giant Branch 

As already mentioned in Sect. [TJ the first evidence that the RGB of NGC6752 is not 
consistent with a single stellar population comes from Grundahl et al. (2002) followed by 
Yong et al. (2008), Milone et al. (2010), Kravtsov et al. (2011), and Carretta et al. (2012). 
These studies have detected a large spread in in the c y Stromgren index^] with the presence 
of three RGBs. The large data set listed in Table 1 allowed us to add further information 
on the multimodal RGB of NGC 6752. 

To extend our multi-wavelength study of the RGB, we attempted to identify the se- 
quences corresponding to each stellar population along the RGB. We already noted that the 
^F275W versus cf275w,F336W,F4iom diagram of Fig. [3] shows a triple RGB. In Fig. [12] panel 
(a) we isolated by hand three groups of stars. Hereafter, we will name these three groups as 
RGBa, RGBb, and RGBc, colored green, magenta, and blue, respectively. The red fiducial 
line is obtained with a procedure similar to the one introduced in Sect. El for MS stars. The 
only difference is that for RGB stars we used a second-order polynomial to interpolate the 
median color and magnitude values measured in the different magnitude intervals. 

The verticalized mF275W versus Acf275W,F336W,F4iom diagram is plotted in panel (b) 



2 The index c y is denned as c y =ci — (b — y), where c\ = (b — v) — (v — b) (Yong et al. 2008). 
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Fig. 11. — Collection of CMDs zoomed around the SGB. SGBa, SGBb, and SGBc stars 
defined in Fig. [TOH are plotted green, magenta, and blue, respectively. 

of Fig. H21 while panel (c) of the same figure shows the histogram of the distribution in 
Acf275w,F336W,F4iom • The histogram is fitted by the sum of three Gaussians which we col- 
ored green, magenta, and blue respectively. From the area of the Gaussians we derive the 
relative fraction of RGBa, RGBb, RGBc stars to be (0.31±0.03, 0.41±0.02, 0.28±0.03). 

In order to get additional information on the three RGBs from all possible combi- 
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F275W, F336W, F410M . 

aC F275W, F336W, F410M 



Fig. 12. — Panel (a): Zoom in of the mp275w versus cf275w,F336W,F4iom diagram of Fig. [3b 
around the RGB. Selected RGBa, RGBb, and RGBc stars are plotted green, magenta, and 
blue respectively. Panel (b): Verticalized diagram for RGB stars with 16.65 < ?tif275w < 
17.95. Panel (c): Histogram of the distribution of stars shown in panel (b). The gray line is 
the best-fitting least-square function defined as the sum of the green, the magenta, and the 
blue Gaussians. 

nations of magnitudes in the photometric passbands of our data set, we follow the same 
procedure for the RGB stars that we performed in Sect. [3] for the MS. Results are illustrated 
in Fig. [T3] where we show the fiducial polynomials in twelve m F si4w versus m x — mpsuw 
CMDs, where X=F225W, F275W, F336W, F390M, F390W, F395N, F410M, F467M, F502N, 
F547M, F555W, and F606W. In this analysis, we did not include F110W and F160W filters 
because, the RGB fiducials are poorly determined in these colors, on account of the small 
number of RGB stars that have IR photometry This is consequence of three main facts: i) 
the IR/WFC3 camera has a smaller field of view than UVIS/WFC3, ii) due to its large pixel 
scale, the IR/WFC3 detector is more affected by crowding, hence the fraction of stars with 
high-accuracy photometry is smaller, iii) IR photometry saturates a couple of magnitudes 
above the MS turn-off. 

Both RGBb and RGBc are typically bluer than the RGBa, with the exception of CMDs 
based on the m F336W — m F si4w an d ?^F39ow — ^F8i4W colors. In the other filters the color 
distance from the RGBa of both the RGBb and the RGBc increases with the color baseline, 
with the possible exception of mF390M— "Ifsuw and mp395N— ^F8i4W- A comparison of Fig. [9] 
and [13] reveals that the behaviour of the three RGBs and the three MSs is quite similar, and 
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Fig. 13. — RGB fiducials in twelve mF8i4w versus mx — mpsuw CMDs (X=F225W, F275W, 
F336W, F390M, F390W, F395N, F410M, F467M, F502N, F547M, F555W, and F606W). The 
color distance from the RGBa, measured at ^pskw = 14.7, is indicated in each panel. The 
positions of RGBa, RGBb, and RGBc at mpg\ 4W =14.7 are represented with green, magenta, 
and blue circles, respectively in the inset of each CMD. 



it will be discussed in Sect. [6j 
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5.1. The chemical composition of the three stellar populations 

In the past three decades many spectroscopic studies have provided us with an accurate 
picture of the chemical composition of NGC6752 (e.g. Norris et al. 1981, Grundahl et al. 
2002, Yong et al. 2003, 2005, 2008, Carretta et al. 2007, 2012). We know that NGC6752 is 
a moderately metal poor ([Fe/H]~ —1.6 Yong et al. 2005, Carretta et al. 2010) GC, with a 
large star-to-star variations in O, N, Na, Mg, and Al. Nitrogen is correlated with aluminum 
and sodium, and has a possible small amplitude correlation with a, Fe-peak, and s-process 
elements (Yong et al. 2008). Both the Na-0 and the Mg-Al anticorrelations have been 
observed by Yong et al. (2005) and Carretta et al. (2007, 2009, 2012). 

Stromgren photometry can provide additional, important information on the chemical 
properties of the stellar populations. Grundahl et al. (2002) and Sbordone et al. (2011) 
have demonstrated that the c y index correlates with the nitrogen abundance and hence can 
be used to identify different stellar populations in GCs. Other combinations of Stromgren 
filters, are sensitive to the chemical differences of GC stars (see, e.g. Grundahl et al. 2000, 
Marino et al. 2011, Carretta et al. 2011). 

The y versus (u — b)—(v — y) diagram of NGC6752 is shown in Fig. [TH upper panel 
(see Carretta et al. 2011 for a discussion on this color index). To avoid the central cluster 
regions, where the photometric error is larger because of crowding, in the upper panel of 
Fig. Ulwe show stars with a distance from the center larger than 1.7 arcmin. As already 
pointed out by Grundahl et al. (2002) in Stromgren photometry CMDs the RGB shows 
three main components. We investigated whether these three RGBs correspond to the three 
populations we have identified in this paper. To do this, we cross-identified the stars in the 
HST and Stromgren catalog and in the CMD in the inset of the upper panel of Fig. [TH (which 
includes only stars out to 1.7 arcmin from the cluster center), we colored green, magenta, 
and blue the RGBa, RGBb, and RGBc stars identified in Fig. [121 Even if the stars common 
to both HST and Stromgren data sets are all located near the cluster center - and hence have 
larger photometric errors - it is clear that stars in RGBa, RGBb, and RGBc have different 
location (color) in the CMD from Stromgren photometry. 

In order to estimate the number of stars in each RGB, we have first isolated the RGB 
region with 13.5 < y < 16.5 and distance from the center larger than 1.7 arcmin, and drew 
by hand a fiducial line through the middle RGB as illustrated in the bottom-left panel of 
Fig. dH Then, we verticalized the RGB following the procedure described in Sect. [3] (middle 
lower panel of Fig. Hlj) . Finally, we obtained the histogram of the A((it — b) — {v — y)) color 
distribution shown in the right panels of Fig. [TH The histogram was least-squared fitted 
with three Gaussians. By summing the area under the three Gaussians, we find that the 
fractions of RGBa, RGBb, and RGBc stars at radial distances between 1.7 and 6.1 arcmin 
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Fig. 14. — Upper panel: y versus (u — v) — (b — y) diagram for NGC6752 obtained using the 
Stromgren photometric catalog from Grundahl et al. (2002). Only stars with distance from 
the cluster center greater than 1.7 arcmin are plotted in this panel. Asymptotic giant branch 
and HB stars, identified in the b versus v — y CMD are represented with gray triangles and 
crosses, respectively. The inset is a zoom around the RGB for stars with distance from the 
center 0.7<R<1.7 arcmin where we colored green, magenta, and blue the RGBa, RGBb, 
and RGBc stars selected in Fig. [12j Lower panels illustrate the procedure to estimate the 
fraction of stars in each RGB. Left panel shows the fiducial line superimposed to the RGBb, 
the rectified CMD for RGB stars is plotted in the middle panel while right panel shows the 
histogram of the A((w — b) — (v — y)) distribution with the best fit Gaussian colored green, 
magenta, and blue. 
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are 27±4%, 44±4%, and 29±3%, respectively. The quoted uncertainties are Poisson errors 
and represent a lower limit of the true errors. All methods applied thus far to the MS and 
RGB, yield very similar fractions of the a, b, and c populations suggesting an association 
between the MS and RGB for the three populations. 

Available spectroscopic abundances allow us to better characterize the chemical content 
of each single RGB. Figure [T5l reproduces the correlations among Al, Mg, N, Na, O using the 
chemical abundance measurements by Yong et al. (2005, 2008). In the upper-right panel, we 
arbitrarily selected three groups of Na-poor, Na-intermediate, and Na-rich stars and colored 
them green, magenta, and blue, respectively. These color codes are used consistently in the 
other panels of this figure. 

As Yong and collaborators have already pointed out, Na-rich stars are enhanced in Al, 
N, and s-process elements, and depleted in O and Mg. Na-intermediate stars are also Al, 
N-enhanced and O-depleted, but their abundance variations are smaller, on average, than 
those of Na-rich stars. There is no significant difference in the Mg and Y content of Na- 
poor and Na-intermediate stars. It is interesting to note that the three groups of stars are 
not chemically homogeneous, as they show star-to-star variation in the abundance of some 
elements that are significantly larger than observational errors (see Yong et al. 2008 for more 
details). 

In the y versus (u — v )-(b — y) diagram of the left-bottom panel of Fig. we mark 
with full dots the Na-poor, Na-intermediate, and Na-rich stars selected in the upper-leftmost 
panel. The three RGBs correspond to groups of stars with different light-elements content, 
as first noticed by Grundahl et al. (2002). We have HST photometry for only one of the 
stars observed by Yong et al. (2008). Its position in the mp275w versus cf275w,F336W,f4iom 
diagram is shown in the lower-middle panel and confirms that the RGBc is made of Na-rich 
stars. 

The average abundance for 22 elements from Yong et al. (2008) for the three groups of 
RGB stars defined above are listed in Table [2j We emphasize how the available spectroscopic 
measures offer us a precious opportunity to characterize the chemical composition of stars 
of the three stellar populations in NGC 6752. 

5.2. Some considerations on the spectroscopic and photometric observational 

evidence of multiple stellar populations 

The Na-0 anticorrelation has been proposed as a possible proof of the presence of 
multiple stellar generations in star clusters. Several authors suggested that the different 
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Fig. 15. — Upper and right panels: Correlations and anticorrelations among the abundances 
of several chemical species from Yong et al. (2003, 2008). We have selected three groups of 
stars with different [Na/Fe] and colored them green, magenta, and blue. In the bottom-left 
and central panel we mark the position of these stars in the v versus (u — v) — (b — y) diagram 
and the mF275W versus cf275w,F336W,F4iom diagram. Stars for which both spectroscopic and 
photometric measurements are available are marked with black circles. 
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Table 2: Average chemical abundance for the three RGB components. 



population 


Population a 


Population b 


Population c 


Element 


abundance [dex] 


a 


N 


abundance [dex] 


a 


N 


abundance [dex] 


a 


N 


[O/Fc] 


0.65±0.04 


0.11 


10 


0.43±0.03 


0.13 


18 


0.03±0.04 


0.11 


10 


[N/Fc] 


-0.11±0.09 


0.25 


8 


0.92±0.13 


0.28 


6 


1.35±0.04 


0.10 


7 


[Na/Fe] 


-0.03±0.02 


0.06 


10 


0.26±0.02 


0.08 


18 


0.61±0.02 


0.05 


10 


[Mg/Fc] 


0.51±0.01 


0.02 


10 


0.49±0.01 


0.02 


18 


0.40±0.02 


0.07 


10 


[Al/Fc] 


0.28±0.05 


0.16 


10 


0.70±0.03 


0.14 


18 


1.14±0.04 


0.12 


10 


[Si/Fc] 


0.27±0.02 


0.07 


10 


0.33±0.01 


0.05 


18 


0.35±0.01 


0.04 


10 


[Ca/Fc] 


0.21±0.03 


0.08 


10 


0.24±0.02 


0.09 


18 


0.27±0.02 


0.06 


10 


[Sc/Fe] 


-0.05±0.01 


0.04 


10 


-0.04±0.01 


0.04 


18 


-0.04±0.01 


0.04 


10 


[Ti/Fe] 


0.10±0.02 


0.06 


10 


0.14±0.01 


0.04 


18 


0.15±0.01 


0.03 


10 


[V/Fc] 


-0.34±0.05 


0.16 


10 


-0.29±0.03 


0.12 


18 


-0.25±0.03 


0.08 


10 


[Mn/Fc] 


-0.50±0.04 


0.11 


10 


-0.44±0.01 


0.06 


18 


-0.45±0.01 


0.04 


10 


[Fe/H] 


-1.65±0.02 


0.07 


10 


-1.61±0.01 


0.02 


18 


-1.61±0.01 


0.01 


10 


[Co/Fc] 


-0.03±0.03 


0.09 


10 


-0.00±0.02 


0.06 


18 


-0.06±0.02 


0.07 


10 


[Ni/Fc] 


-0.06±0.02 


0.07 


10 


-0.06±0.01 


0.04 


18 


-0.03±0.01 


0.03 


10 


[Cu/Fc] 


-0.66±0.03 


0.09 


10 


-0.59±0.01 


0.05 


18 


-0.60±0.01 


0.04 


10 


[Y/Fc] 


-0.09±0.02 


0.06 


10 


-0.01±0.02 


0.06 


18 


0.01±0.01 


0.02 


10 


[Zr/Fc] 


0.07±0.05 


0.15 


10 


0.20±0.02 


0.07 


18 


0.21±0.03 


0.09 


10 


[Ba/Fe] 


-0.09±0.04 


0.13 


10 


-0.12±0.03 


0.13 


16 


0.05±0.02 


0.07 


9 


[La/Fc] 


0.12±0.02 


0.02 


2 


0.10±0.01 


0.04 


12 


0.13±0.04 


0.06 


3 


[Ce/Fe] 


0.28±0.03 


0.09 


10 


0.25±0.01 


0.04 


18 


0.28±0.02 


0.06 


10 


[Nd/Fe] 


0.23±0.01 


0.04 


10 


0.22±0.01 


0.05 


18 


0.23±0.01 


0.04 


10 


[Eu/Fe] 


0.31±0.03 


0.10 


10 


0.30±0.02 


0.08 


18 


0.34±0.03 


0.10 


10 
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populations of stars in GCs can be identified on the basis of their position along the Na-0 
abundance plane, with the first generation including stars with oxygen and sodium abun- 
dance similar to halo field stars at the same metallicity, with remaining Na-rich/O-poor to 
be considered as part of a second generation (e.g. Kraft et al. 1994). 

Spectroscopic studies show that the Na-0 anticorrelation is a common property among 
GCs (e.g. Ramirez & Cohen 2002, Carretta et al. 2009, 2010 and references therein). How- 
ever, there clearly is a problem we have not solved yet. Figure 12 from the compilation by 
Ramirez & Cohen (2002) and Fig. 1 by Carretta et al. (2010) shows an almost continuous 
distribution of stars in the Na versus O plane, despite the fact that photometric evidence 
in many of the clusters included in that figure (e.g. NGC 2808, NGC 6397, NGC 6752 as 
shown in the present paper) show multimodal, maybe discrete, sequences in the CMD when 
high-accuracy photometry on images collected with the appropriate filters is used^. Why do 
we have this possible difference between the spectroscopic and the photometric manifestation 
of the multiple stellar populations in GCs? 

Carretta et al. (2009) suggested criteria to separate stellar populations on the basis 
of their position in the Na-0 plane. They defined as primordial (P) component all stars 
with [Na/Fe] ratio in the range between [Na/Fe] m i n and [Na/Fe] m i n +0.3 where [Na/Fe] m i n 
is the minimum value of the ratio [Na/Fe] ratio estimated by eye. The remaining stars 
are considered all second-population stars, and have been further divided into two groups. 
Stars with ratio [0/Na]> —0.9 dex belong to the intermediate (I) population, while those 
with [0/Na]< —0.9 dex are defined as extreme (E) population. Clearly, this separation is 
arbitrary and not based on any feature (gaps or peak) in the NaO diagram, and it has no 
clear physical meaning. As an example, in Fig. [161 we a Pply the same criteria as proposed 
by Carretta et al. (2009) to the sample of NGC 6752 stars studied by Yong at al. (2008) and 
further analyzed in Sect. 15.11 The two red segments define the regions in the Na-0 plane 
populated by P, I, and E stars, and are determined following the recipes by Carretta and 
collaborators. The stars of the populations a, b, and c identified in this paper are colored 
green, magenta, and blue, respectively. Figure [TBI shows that: i) the Carretta et al. (2009) P 
component includes all the stars in population a but is contaminated by population b stars; 
ii) the I component contains both population b and population c stars, and (iii) no stars 
belong to the E component. Apparently, the general criteria introduced by Carretta et al. 
(2009) do not apply to NGC 6752. 



3 Note however that, in some cases, multimodal distribution in Na and O have been detected also from 
high-resolution spectroscopic (see e.g. Yong et al. 2008, Marino et al. 2008, Lind et al. 2011 for the cases of 
NGC 6752, M4, and NGC 6397). 
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The main problem here is not related to the meaning of the Carretta et al. (2009) 
definition. The problem is that there still is an inconsistency between what we observe with 
spectroscopy (a continuous distribution along the Na vs. O plane) and what accurate, high 
precision photometry tells us, i.e. that most clusters host distinct, separate evolutionary 
branches in the CMD. It is important to understand whether this is just the consequence of 
the (internal) errors in the measurement of Na and O abundances, or whether there is some 
underlying physical reason we have not understood, yet. 




Fig. 16. — Reproduction of the Na-0 anticorrelation from Yong et al. (2008). Stars belonging 
to the populations a, b, and c defined in this paper are colored green, magenta, and blue, 
respectively. The red segments indicate the separation for the P, I, and E components 
suggested by Carretta et al. (2009). 



6. The Helium abundance of the three stellar populations of NGC 6752 

In this section, we will use the multicolor HST photometry to further characterize 
the chemistry of the different stellar populations in NGC 6752 and estimate the helium 
abundances of its three stellar populations. The multi-dimensional space of our CMDs 
makes it difficult to visualize fitting isochrones to our sequences. To make the comparison 
easier, we will quantify the color separation between the sequences at two fiducial points at 
the level of the MS and the RGB. 

In all fourteen CMDs of Fig. [91 we calculated the color differences between MSa, MSb, 
and MSc at the reference magnitude mpg\ 4W =18.5. These color differences are calculated by 
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subtracting from the color of the MSa fiducial at mpg 14W =18.5 the color of the MSb (or MSc) 
fiducial at the same luminosity. Left panel of Fig. [T7] shows the measured color difference as 
a function of the central wavelength of the mx filter. 




5000 10000 15000 2000 4000 6000 8000 

central X[A] central X[A] 

Fig. 17. — Left Panel: mx—mFsuw (or mpsuw-wx) color distance between MSb and MSa 
(magenta circles) and between MSc and MSa (blue triangles) as a function of the central 
wavelength of the X filter. Right Panel: Color distance between RGBb and RGBa (magenta 
circles) and between RGBc and the RGBa (blue triangles). The color distances of MS and 
RGB sequences are measured at the reference magnitude mpg 14W =18.5 and m F g 14W =14.7, 
respectively. 

It is now clear why the cf275w,F336W,F4iom index is an efficient tool to identify multiple 
sequences in the CMD. Figure [T7] indeed shows that both the mp275w — m F336W an d the 
m F336W — m F4ioM color provide large separations between the MSa (or the RGBa) and the 
other two MSs (RGBs). The MSa (RGBa) is redder than the MSb and the MSc (RGBb 
and RGBc) in m F 275w — m F336W, but it moves to the blue in mp336w — m F4ioM- Therefore, 
the color index c F2 75W,F336W,F4ioM = (mF275W - ^F336w) - ("1F336W - ™F4iom) maximizes the 
separation among the sequenceso 



4 By following a similar approach we suggest that other color indices cf336W,F4iom,F814W=(wf336W — 
wf4iom) — (wf4iom — WF8i4w) can be powerful tools to detect multiple stellar populations in the CMD of 
GCs. The later is less efficient than cf275W,F336W,F4iom, but it uses a wavelength range that is accessible to 
ground-based telescopes. 



-28 - 



In order to compare these observations with expectations from synthetic photome- 
try, we followed a procedure that has already been used in previous publications (Milone 
et al. 2012a,b, Bellini et al. 2013, see also Sbordone et al. 2011). Briefly, we used the 
BaSTI isochrones (Pietrinferni et al. 2004, 2009) for the populations listed in Table [31 
and determined T e ff and logg for a MS star with m^ 14W =18.5. To do this, we assumed 
EiB — V) = 0.03 and (m — M)v = 13.19 in agreement with the values of reddening and dis- 
tance modulus listed in the Harris (1996, updated as in December, 2010). For ACS/WFC fil- 
ters we used the extinction coefficients tabulated by Bedin et al. (2005) for a cold (T=4000K) 
star. For the WFC3/UVIS filters, we linearly interpolated the coefficients tabulated by Bedin 
et al. (2005), and for WFC3/IR we adopted the values listed in the York Extinction Solver 
website^. For each of the three populations, we adopted the abundances of the 22 elements 
listed in Table [2J which represent the average abundance of the three stellar populations. 
Since carbon abundance was not measured by Yong et al. (2008) we adopted the [C/Fe] 
values from Carretta et al. (2005) and take [C/Fe]=0.15 for population a, and [C/Fe]=— 0.15 
for both population b, and c. 

We assumed that the population a has primordial helium content (Y=0.246), and 
adopted for both the population b, and c different helium abundances as described in the 
following. We used the ATLAS 12 code (Kurucz 2005, Castelli 2005, Sbordone 2007) to 
calculate model atmospheres by using the specific chemical composition of each stellar pop- 
ulation, and assuming the temperatures and gravities listed in Table |3j Then, we used the 
SYNTHE code (Sbordone et al. 2007) to synthesize the spectrum from 1000A to 20 000A 
and the resulting spectra were convolved with the transmission curves of the system formed 
by the telescope, the camera, and each of our filters to produce the synthetic magnitudes 
and colors for each photometric band. 

The adopted chemical compositions for the three different stellar populations a,b,c iden- 
tified within NGC6752 , (MSs, SGBs, and RGBs), are tabulated in Table[3]for three different 
options. In all these three options, population a is assumed to have a canonical helium abun- 
dance and the chemical composition in Table 2. Instead, the chemical compositions and the 
helium fraction (Y) of populations b and c, are different for different options. In details: 

• In Option I, the three populations have the same chemical composition, but three 
different values for the helium content (Y=0.246, 0.254, and 0.275). 

• In Option II, the three populations have the same helium fraction, but the three chem- 
ical compositions given in Table 2. 



5 http://www3.cadc-ccda.hia-iha.nrc-cnr.gc.ca/community 
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• In Option III, the three populations have the three chemical composition given in 
Table 2, and three different values of helium as described below. We estimated for 
population a, b, and c, Y=0.246, 0.254, and 0.275, respectively. 

The comparison of the observed mx — mpsuw differences between MSa and MSc against the 
synthetic ones are shown in the left panel of Fig. [TH 

The blue squares indicate the color differences corresponding to Option I, where we 
assumed for the two stellar populations the same element abundance as for MSa, but a 
different helium content. Similar to what observed in the cases of 47Tuc and NGC6397, we 
find a good agreement in most bands but a significant disagreement for filters bluer than 
~4000A, and conclude that helium cannot be the only parameter responsible for the color 
differences between MSa and MSc. This confirms that the observed light-element abundance 
of the three stellar populations play a fundamental role in the MS morphology. 

In Option II we assumed for MSc stars the same He content as for the MSa, and the 
chemical composition listed in Table [31 The colors resulting from this option are represented 
by gray triangles in Fig. [TSJ and indicate a significant disagreement between observed and 
synthetic colors. It is worth noting that the observed light-element difference among the 
three stellar populations have a negligible effect for filters redder than ~4000Aas already 
shown by Sbordone et al. (2011). 

Finally, in Option III (red asterisks) we assumed differences in both helium and chem- 
ical abundances. To determine the value of Y that best reproduces the observed points, 
we generated a grid of synthetic spectra by assuming different helium abundances, with Y 
ranging from 0.246 to 0.290 in steps of AY=0.001. For each synthetic spectrum we deter- 
mined the mx — mpsuw color distance between MSc and MSa corresponding to each X filter: 
(Acolor syn ), and calculated d{Y) = \ Acolor syn — Acolor| where Acolor is the observed color 
distance between the two MSs in that filter. 

The helium abundance (Y*) that minimizes d(Y) is assumed as the best estimate of 
Y for the filter X. The helium content of MSc is then calculated as the weighted mean 
of the available Y* measurements by using only those filters redder than F395N we obtain 
Y=0.273±0.002. The quoted error comes from the weighted mean, and does not take into 
account neither the uncertainties of the synthetic spectra, nor possible errors on the value 
we assumed for the primordial helium. Note that we have excluded from this analysis 
all the UV and far-UV filters as they are very sensitive to small variations of light-element 
abundance. Option III (red asterisks), provides the best agreement with observations though 
not completely satisfactory at short wavelengths. 

The right panel shows the comparison of the observed nix ~ m -F8i4W color differences 
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between the MSa and the MSb against the synthetic one. Again in this case the best fit is 
given by the Option III, and the agreement between synthetic and observed photometry is 
much better than what we obtained for the MSc and MSa color differences. By using the 
procedure described above we obtain for MSb Y=0.253±0.001. 

The synthetic spectra of a MSa and a MSc star, as calculated for the Option III, are 
shown in the upper panel of Fig. [T9j The difference between the two spectra is in the middle 
panel, while the lower panel shows the (normalized-to-peak) band-passes of the filters used 
in this paper. 

In summary, the observed color differences between the three MSs are consistent with 
three populations with different helium and light-element abundances. Specifically the MSa 
corresponds to the first stellar population with primordial He and O-rich/N-poor stars, 
the MSc is made of stars enhanced in He and N, but depleted in O, and MSb stars has 
intermediate He and C, N, O abundances. 
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Fig. 18. — Color separations at different color baselines for the ridge lines of the MSc and 
the MSa (left), and the MSb and the MSa (right) at mpsuw =18.5. Observations are plotted 
as blue triangles and magenta circles for the MSc and the MSb, respectively, while the color 
differences expected from theoretical Options I, II, and III are shown as blue squares, gray 
triangles, and red asterisks. On the right of each panel, the small boxes correspond to regions 
centered on the F275W, F336W, F410M, and F555W bands. 

In Fig. [201 for the three RGBs we repeated the same analysis as done for the MSs. We 
have calculated the color distance of RGBb and RGBc from RGBa at WmVw =14.7. In 
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Fig. 19. — Upper Panel: comparison of the synthetic spectra of a MSa star (green) and a 
MSc stars (blue). Middle Panel: Difference between the spectra of the MSa and the MSc 
star. The location of relevant molecular features are indicated. Lower Panel: Normalized- 
to-peak response of the HST filters used in this paper. 
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Fig. [201 these color differences are plotted as a function of each filter and labeled with its 
corresponding effective wavelength. 

Observations are compared with synthetic colors calculated for a RGB stars with rap8 14W =14.7 
by using the same procedure described for MS stars, and for the same three Options I, II, 
and III listed in Table [3j Similar to what found for the MSs, the color differences between 
the three RGBs can be best reproduced by Option III, which reinforces the idea previously 
proposed that RGBa, RGBb, and RGBc stars are the progeny of MSa, MSb, and MSc stars, 
respectively. Specifically, by using the procedure described above for the MSs we obtain for 
RGBc and RGBb helium abundances of Y=0.272±0.005, and Y=0.255±0.004, respectively. 




Fig. 20. — Left Panel: mx- ^fsuw color separation between RGBc and RGBa as a function 
of the central wavelength of the corresponding X filter. Right Panel: Color separation 
between RGBb and RGBa. The color distances are measured at the reference magnitude 
m F8i4W~14-7- Observations are plotted as blue triangles for RGc and magenta circle for 
RGBa, while the color differences expected from theoretical Options I, II, and III are shown 
as blue squares, gray triangles, and red asterisks. The small boxes correspond to regions 
centered on the F275W, F336W, F410M, and F555W bands. 



-33 - 



Table 3: Parameters used to simulate synthetic spectra of an MSa, MSb, and an MSc star 
with mF8i4W=18.5, and an RGB star with mF8i4W=14.7 for the three assumed options. For 
all the populations we assumed [Fe/H] = —1.6. The adopted chemical composition is given 
in Table |2J 



MS (Option) 


T c ff 


logs 


Y 


chemical composition (see Table [2]l 


MSa (all) 


5445 


4.68 


0.246 


as for population a 


MSb (I) 


5470 


4.68 


0.254 


as for population a 


MSb (II) 


5445 


4.68 


0.246 


as for population b 


MSb (III) 


5470 


4.68 


0.254 


as for population b 


MSc (I) 


5534 


4.69 


0.275 


as for population a 


MSc (II) 


5445 


4.68 


0.246 


as for population c 


MSc (III) 


5534 


4.69 


0.275 


as for population c 


RGBa (all) 


5343 


3.26 


0.246 


as for population a 


RGBb (I) 


5351 


3.26 


0.254 


as for population a 


RGBb (II) 


5343 


3.26 


0.246 


as for population b 


RGBb (III) 


5351 


3.26 


0.254 


as for population b 


RGBc (I) 


5373 


3.25 


0.275 


as for population a 


RGBc (II) 


5343 


3.26 


0.246 


as for population c 


RGBc (III) 


5373 


3.25 


0.275 


as for population c 



6.1. The radial distribution of the three stellar populations 

In order to calculate the radial distribution of the three stellar populations in NGC 6752, 
we first divided our sample into two regions. The first one includes stars with radial distance 
from the cluster center smaller than 1.7 arcmin, and is covered by HST observations. The 
second one contains stars with radial distances from the center between 1.7 and 6.5 arcmin, 
which are included in the ground-based catalog. The upper panel of Fig. [21] shows the 
fractions of MSa, MSb, and MSc stars with respect to the total number of MS stars as green, 
magenta, and blue triangles, and the fractions of RGBa, RGBb, and RGBc stars with respect 
to all RGB population as green, magenta, and blue circles. 

Then, we further divided the region with radial distance smaller than 1.7 arcmin into 
four circular subregions, each one containing almost the same number of MS stars. For 
each region, we calculated the fraction of MSa, MSb, and MSc stars by using the procedure 
already described in Sect. |3j Similarly, we divided the region with radial distance between 
1.7 and 6.5 arcmin into two parts, containing almost the same number of RGB stars and 
calculated the fraction of RGBa, RGBb, and RGBc stars as in Sect. 15.11 Results are shown 



6 We note here that the fraction of stars in the different MSs has been estimated by comparing stars in 
small magnitude intervals (see Sect. 3). Since these stars have similar luminosities and the color differences 
between MSa, MSb, and MSc stars is quite small, incompleteness is not an issue in this comparison. In any 
case, the completeness level is greater than 0.67 for the adopted magnitude interval, with mpgi4vv < 20.3. 
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in the lower panel of Fig. |2TJ Note that, due to the relatively small number of RGB stars in 
the HST field we preferred to limit the analysis based on HST data to MS stars. 

The fraction of population a, population b, and population c stars are listed in Table 4, 
where we also included the minimum and the maximum radius of the circular region (R m m 
and -Rmax), and the average radial distance of the MS or RGB stars used to estimate the 
populations ratio (R ave ). The later is calculated as the mean radius of the MS or RGB stars 
in that bin. 

There is no significant radial trend in the relative numbers of the three stellar populations 
within six arcmin from the cluster center. Apparently these results do not support the 
conclusions by Kravtsov et al. (2011) who, analysed ground-based images, and reported a 
strong difference in the radial distribution between the RGB subpopulations. These authors 
found that the change in the fraction in the two RGBs occurs at radial distance close to 
the half-mass radius of the cluster (rh=2.34 arcmin, Harris 1996, 2003) and becomes much 
stronger at larger radial radii. An extension of the analysis presented in this paper to larger 
radial distances from the cluster center is mandatory to properly characterize the radial 
distributions of the different subpopulations. 

As pointed out by the referee, the relaxation time of NGC 6752 at the time the secondary 
generations formed is a fundamental ingredient for a proper interpretation of these results. 
The estimate of this quantity is beyond the objectives of this paper and is complicated 
by the fact that NGC 6752 (as, in general, all GCs showing evidence of multiple stellar 
populations) may have been significantly more massive at the time of the formation of their 
stellar generations, (see e.g. Conroy et al. 2012, Goudfrooij et al. 2011, D'Ercole et al. 2010). 
The populations ratio listed in Table 4 at several radial distance from the cluster center can 
provide useful contraints for the models of formation and evolution of stellar populations in 
GCs. 



7. Summary 

We analyzed HST images and ground-based catalogs obtained through a large set of 
filters in order to identify multiple stellar populations in NGC 6752. We find that the MS 
of this cluster splits into three components, in close analogy with what we observed for its 
RGB and its SGB. We conclude that NGC 6752 hosts at least three stellar populations, 
whose evolution can now be followed from the MS up to the RGB tip. This result is nicely 
summarized in Fig. [221 where we show some representative CMDs where each population 
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Fig. 21. — Radial distribution of the fraction of population a (green symbols), population 
b (magenta symbols), and population c (blue symbols) with respect to the total number of 
(population a + population b + population c) of stars. Circles and triangles refer to the 
measures obtained from RGB and MS stars, respectively. The dotted and the dashed-dotted 
vertical lines mark the core and the half-mass radius, respectively. In the upper panel we 
have used one single radial interval for HST and one for ground-based data. In the lower 
panel we have divided the HST field of view into four radial bins, and the region with radial 
distance larger than 1.7 arcmin (ground-based data) into two bins. See text for details. 
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-^min 


Rmax 


Rave 




population ratio 




sequence 








a 


b 


c 




0.00 


1.70 


0.95 


0.25±0.02 


0.44±0.04 


0.31±0.03 


MS 


0.00 


1.70 


0.87 


0.28±0.03 


0.41±0.03 


0.31±0.03 


RGB 


1.70 


6.13 


3.26 


0.27±0.04 


0.44±0.04 


0.29±0.04 


RGB 


0.00 


0.53 


0.31 


0.24±0.02 


0.47±0.05 


0.29±0.042 


MS 


0.53 


0.83 


0.68 


0.23±0.02 


0.52±0.04 


0.25±0.031 


MS 


0.83 


1.12 


0.97 


0.28±0.02 


0.42±0.05 


0.29±0.042 


MS 


1.12 


2.33 


1.44 


0.28±0.03 


0.44±0.04 


0.28±0.025 


MS 


1.70 


3.11 


2.35 


0.26±0.05 


0.45±0.05 


0.29±0.05 


RGB 


3.11 


6.13 


4.15 


0.30±0.05 


0.43±0.05 


0.27±0.05 


RGB 



Table 4: Fraction of population a, population b, and population c stars calculated in different 
circular regions with different radial distance from the cluster center. The minimum and 
maximum radius (-R m i n and -R max ) of each region are listed together with the average radial 
distance of stars in each region (i? avB ). The last column indicates the sequence of the CMD 
(MS or RGB) used to estimate the populations ratio. 

can be followed along its evolutionary phasesJll 

The multi-wavelength photometric data set allowed us to complement and extend the 
information on the chemical composition of the different populations available from spectro- 
scopic chemical abundance measurements of a limited sample of stars. We calculated model 
atmospheres for MS and RGB stars accounting for the available chemical composition, and 
demonstrated that the three groups of stars have different helium and light-element abun- 
dances. The most straightforward interpretation is that Group 'a', which contains about 
25% of the total number of stars, is the first stellar population, originated in a molecular 
cloud with a chemical composition similar to that of the Galactic halo, of which it shares its 
chemical composition. The majority of stars represents a second (and third) stellar popula- 
tions which we named 'b' (and 'c') and contain about 30% (and 45 %) of the cluster stars. 
They formed out of material that had been partly processed through first-generation stars, 
and are C/O poor, N/Na/Al rich, and enhanced in helium by Ay ~0.01 and ~0.03, respec- 
tively. We note that stars in each RGB sequence exhibit a wide spread in the abundance of 
some light-elements (e.g. N and Mg), thus suggesting that the three groups of stars defined 
above are not chemically homogeneous. Interestingly also the HB morphology of NGC 6752 
seems to be composed also by three sub-groups (Momany et al. 2002), suggesting that these 
three populations eventually evolve in HB stars with different properties. No radial gradient 



The fiducial lines of the three stellar populations are available at this url: 
http: / / www.astro.unipd.it /globulars / 
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Fig. 22.— m F27 5W versus m F275W -m F 33 6W (top) and m F336W versus m F336W -m F 395 N CMDs. 
We have colored green, magenta, and blue the three groups of stars selected in Figs. (5J [101 
and [121 The fiducial lines of the three stellar populations of NGC 6752 are shown on the 
right. 
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of the different stellar populations was detected within NGC6752. 

Similar conclusions have been reached for other GCs: the multiple MSs of u Centauri and 
NGC2808 suggest extreme helium abundances (Y ~ 0.39, e.g. Bedin et al. 2004, D'Antona et 
al. 2005, Piotto et al. 2007, King et al. 2012) while the multiple sequences of NGC 6397 and 47 
Tuc imply small helium enhancements (AY ~ 0.01—0.02, Di Criscienzo et al. 2010a,b, Milone 
et al. 2012a,b). Interestingly, while in NGC 6397 and 47 Tuc there is evidence of two groups 
of stars with slight different helium content in the case of both NGC 6752 and NGC 2808 
we have identified at least three stellar populations. It is worth noting that a multimodal 
RGB with at least three components has been also observed in NGC 6205 (Grundahl et al. 
1998). A multi band study of this cluster is mandatory to estimate the helium content of its 
stellar populations. Our results provide evidence that differences in helium abundance are 
a quite common feature of different stellar populations in GCs. D'Ercole et al. (2010, 2012) 
show that different combinations of helium and CNO variations can be obtained, due to the 
different modalities in which in each cluster the ejecta of stars are diluted with primordial 
gas. Further study of color-magnitude diagrams multiple populations may allow to constrain 
for each cluster the parameters necessary to its specific model. 

Finally, we discuss a potential disagreement between spectroscopic and photometric 
evidence of multiple stellar populations in GCs which needs to be solved for further progress 
in our understanding of GC stellar populations. 
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